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Abstract 

We point out novel consequences of general relativity involving tidal 
dynamics of ultrarelativistic relative motion. Specifically, we use the 
generalized Jacobi equation and its extension to study the force-free 
dynamics of relativistic flows near a massive rotating source. We show 
that along the rotation axis of the gravitational source, relativistic 
tidal effects strongly decelerate an initially ultrarelativistic flow with 
respect to the ambient medium, contrary to Newtonian expectations. 
Moreover, an initially ultrarelativistic flow perpendicular to the axis 
of rotation is strongly accelerated by the relativistic tidal forces. The 
astrophysical implications of these results for jets and ultrahigh energy 
cosmic rays are briefly mentioned. 

PACS numbers: 04.20. Cv, 97.60.Lf, 98.58.Fd, 98.70. Sa 
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1 Introduction 

This paper is about the motion of free test particles in the background grav- 
itational field of a Kerr black hole. The geodesies of Kerr spacetime are well 
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known and have been extensively studied in standard systems of coordi- 
nates such as the Boyer-Lindquist coordinate system. The Boyer-Lindquist 
coordinates are adapted to the stationary and axisymmetric nature of the 
source and are Minkowskian (expressed in spherical polar coordinates) in- 
finitely far from the black hole. We are interested however in the relative 
motion of free nearby particles; as explained in detail in [T' , from the stand- 
point of relativity theory the study of relative motion is certainly preferred, 
since it is in keeping with the spirit of the theory. Moreover, relative motion 
has direct observational significance. Thus theoretical results obtained in 
the study of relative motion can be directly compared with suitable experi- 
mental data ^El- 

To study the relative motion of free particles, we need a reference particle 
whose motion is well known. For the purposes of the present work, we 
select the free motion of particles on radial escape trajectories along the 
rotation axis of the Kerr black hole. Thus at a given initial position along 
a radial direction in Schwarzschild spacetime or the rotation axis in Kerr 
spacetime, the reference particle has an escape velocity such that the particle 
eventually reaches infinity with zero kinetic energy. We then study the 
relative motion of the relativistic particles that start from the same location 
as the reference particle and travel in different directions. The origin of 
these relativistic particles on the axis of rotation will not be discussed in 
this paper. We will simply assume that such particles exist near the poles 
of the Kerr black hole. The nature of the central engine and the details of 
the accretion phenomena that could generate such relativistic particles are 
beyond the scope of this work. Instead, we concentrate on the gravitational 
dynamics of these particles relative to the reference particle. Moreover, it 
turns out that the existence of event horizons has no bearing on our main 
results; therefore, our treatment applies equally well to test particles in the 
exterior of a gravitational source whose field is adequately described by the 
exterior Kerr spacetime. However, our results are in general observationally 
significant only in the case of highly collapsed gravitational sources such as 
neutron stars and black holes. 

To characterize the relative motion, we establish a quasi-inertial coordi- 
nate system along the worldline of the reference particle. This Fermi normal 
coordinate system [3. is valid in a cylindrical region of radius TZ along the 
worldline of the reference particle. Here TZ is the radius of curvature of the 
spacetime manifold. The geodesic equations of motion for the free particles 
are then integrated in this Fermi coordinate system subject to the limita- 
tion that the relative distance must remain well within the admissible range 
of Fermi coordinates. The results of such an investigation of the relative 
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motion of free particles would naturally depend on the choice of the refer- 
ence particle. For instance, a separate investigation would be necessary if 
the reference trajectory is an orbit in the equatorial plane of the Kerr black 
hole. In fact, a judicious choice of the reference particle would be essential 
to explain a particular observational result. Moreover the relative distance 
in our approach should be less than TZ. 

We find that ultrarelativistic flows in the gravitational field of collapsed 
configurations can exhibit purely non-Newtonian behavior due to general 
relativistic tidal effects. These include the phenomenon of the tidal decel- 
eration along the axis of rotation and the phenomenon of tidal acceleration 
perpendicular to the axis of rotation, which follow from the integration of the 
generalized Jacobi equation We study this equation and its extension 

in this paper and compare our results with the corresponding nonrelativistic 
tidal effects. In the following we use units such that c = 1. The signature 
of the spacetime metric is assumed to be +2. 

Imagine a geodesic worldline in the gravitational fleld of an external 
source. Let w'^ = dx^ /dr be the four- velocity vector of this reference geodesic 
and T be its proper time. A neighboring geodesic with an arbitrary velocity 
relative to the reference geodesic is connected to it by the deviation vector 
^^{t) such that £,^u^ = 0. The generalized Jacobi equation |T| is given by 
the geodesic deviation equation to first order in namely, 



+ 2R%,^-^Cu" + i;R%ua^e^ = 0. (1) 



Dt^ +3 P''' Dt^ Dt 
Here D^i^/Dt = i,^. is the covariant derivative of along the refer- 
ence worldline. Let A'^j.^^ be an orthonormal tetrad frame that is parallel 
propagated along the reference trajectory such that A'^^q^ = u^. Then, 
^P = X^\^^^ and equation (HJ may be written in terms of Fermi coordinates 
(T, X), where T = t along the reference geodesic, as jHE] 



d^X 

- + " RoiojX-' - 

2 rr„ « 2 



^ + ^RoiojX^ + 2 ""RikjoV'X^ 



+(2 ^RokjoV'V' + - ^RikjiV'V' + - ^RokjeV'V^V')X^ = 0. (2) 
Here = dX''/dT, |V| < 1 along the reference geodesic and 

^RalS'jsiT) = R^uptj^^ (a)'^^ iP)'^'' {-y)'^'' {&) ('^) 
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is the projection of the Riemann curvature tensor on the orthonormal tetrad 
of the reference observer along its worldhne. If the relative velocity is much 
smaller than the speed of light, |V| <C 1, the velocity-dependent terms in 

and ((21) may be neglected and the generalized Jacobi equation reduces 
to the standard Jacobi equation that expresses the linear evolution of the 
deviation between two neighboring geodesies that have negligible relative 
motion. This relative motion generally evolves slowly on the scale of the 
radius of curvature TZ, which is defined such that 7^~^ is the supremum 
of l^Ra/s-ysl', therefore, |X| <C TZ in Q, since higher-order terms in the 
deviation equation have been neglected. In the generalized Jacobi equation, 
the relative velocity could approach the speed of light; therefore, the relative 
motion could evolve rapidly thus further restricting the temporal domain of 
validity of equation (jSJ. 

Equation can be derived in a straightforward manner using the re- 
duced geodesic equation in Fermi coordinates. That is, starting from the 
geodesic equation for a free particle in arbitrary coordinates = 

where ds'^ = g^ydx^dx'^ , one can simply derive the reduced geodesic equa- 
tion [Tl 

d'^x^ ~ Q dx"^ dx^ dx^ - ^ dx" dx^^ 

IP " ^^"^^^^^ + ^"^^IT = °- 

Specializing this equation now to the Fermi coordinates established along 
the reference trajectory, i.e. {t,x^) i— > (T, X*), and using 

goo = -l-'"RoiOj{T)X'X^ + ■■■ , 

goi = -'^^Rojik(.T)X^X^ + -- - , 

9ij = 5,, - ^ ^Rikji{T) X'^X' + • • • , (6) 

one arrives at the tidal equation for the relative motion, since in these coordi- 
nates the reference particle remains fixed at X = 0. Restricting our attention 
to the terms given explicitly in © and thus neglecting higher-order terms of 
the metric in the spatial Fermi coordinates, the tidal equation reduces to the 
generalized Jacobi equation We study the main physical consequences 
of this equation in the field of a Kerr black hole (see also Appendix A ). 



Higher-order tidal accelerations are discussed in Appendix B 



Tidal deceleration along the rotation axis is discussed in the next sec- 
tion. The non-Newtonian character of our results is emphasized in section [HI 
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The analysis of the generahzed Jacobi equation is extended to the three- 
dimensional case in section^ Tidal acceleration normal to the rotation axis 
is discussed in section |H1 Finally, section IHl contains a brief discussion of our 
results. 

2 Tidal deceleration 

Imagine a Kerr source of mass M and angular momentum J = Maz at 
the origin of an asymptotically inertial coordinate system {t,x,y,z). Here 
a > is the specific angular momentum of the source. The metric of Kerr 
spacetime can be written in Boyer-Lindquist coordinates as 

ds"^ = -dt^ + S(^dr2 + _^ (^2 _^ ^2^) g-j^2 ^ ^^2 

+ 2GM^{dt-aW{}d<i)f, (7) 

where S = r'^+a? cos^ i? and A = r"^ — 2GMr+a? . Asymptotically (r oo), 
the spherical polar coordinates (r, 0) are related to the Cartesian coordi- 
nates (x, y, z) in the standard manner. In this paper, we are interested in 
the motion of free test particles in the stationary and axisymmetric exterior 
Kerr spacetime. Let us first consider geodesic motion along the axis of rota- 
tion of the source, i.e. the z-axis. In terms of the Boyer-Lindquist temporal 
coordinate t and radial coordinate r, the motion of a free test particle along 
the rotation axis is given by 

dt _ + dr _ 2GMr 

'd^~^r'^-2GMr + a?' d^~^\^ r^ + a^' 

where r is the proper time along the path and 7 > is a constant of 
integration. For 7 < 1, the test particle cannot escape the gravitational 
attraction of the source. When the test particle can escape to infinity, 7 > 1 
and 7 is its Lorentz factor at spatial infinity. In the limiting case of 7 = 1, the 
test particle approaches infinity with zero speed. We assume in the following 
that there are many such outward moving test particles with 7 = 1 in the 
neighborhood of the rotation axis of the Kerr source forming an ambient 
medium around the rotating source. Relative to this ambient medium, we 
wish to study the motion of an ultrarelativistic jet emitted outward from a 
region near the source along the axis of rotation. 

Consider a jet clump starting from tq and moving rapidly along the z- 
axis according to equation Q with 7 > 1. This relative motion is properly 
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described in the first approximation using the generalized Jacobi equation 
expressed in Boyer-Lindquist coordinates. To compare the consequences 
of such an equation with astronomical data regarding the motion of jets, it 
is appropriate to transform this equation to a quasi-inertial Fermi normal 
coordinate system that is established along the reference trajectory. The 
explicit relationship between the Boyer-Lindquist coordinate system of Kerr 
spacetime and the Fermi system is not needed for the derivation of the 
equation of relative motion in Fermi coordinates ((21). We therefore proceed 
keeping in mind that the radial coordinate r of the jet clump is a function of 
the Fermi coordinates, r = r{T,Z), but that this relationship shall remain 
implicit throughout this work. 

Along the worldline of a generic outward moving test particle with 7 = 1 
on the rotation axis, we establish a Fermi coordinate system {T,X,Y, Z) 
based on an orthonormal nonrotating (i.e. Fermi- Walker transported) tetrad 
frame. The free test particle follows a geodesic; therefore, Fermi-Walker 
transport reduces to parallel transport in this case. The Fermi coordinates 
of the reference particle are then (r, 0, 0, 0), i.e. the reference particle is at 
the spatial origin of the Fermi coordinates. The quasi-inertial Fermi system 
is the one appropriate for the comparison of the theory with observation as 
well as Newtonian physics [2- In fact, the speed of a jet is determined by 
monitoring the motion of a jet clump with respect to neighboring "fixed" 
features of the ambient medium [3]. The motion of the jet clump relative to 
the reference particle in the Fermi system is given by the tidal equation 

f^ + k{l-2V')Z + OiZ')=0, (9) 

where V = dZ/dT and the curvature k = '^iioaos is given by 

^r(r2 -3a2) 
(r^ -|- 

Here r is the Boyer-Lindquist radial coordinate of the reference particle and 
r(T) can be determined from the integration of the second equation in Q 
with 7 = 1, r = r and dr/dT > 0. We will be concerned with regions along 
the rotation axis such that 3M^ > 3a^, so that k < 0. Solving equation 
(jH)) with 7 = 1 and t ^ T for r(T) and substituting this in equation (|lfl|) 
gives us the appropriate k(T) to use in the tidal equation ©. Let us note 
that equations (|g|)-(|lfl|) are valid for any free reference particle moving along 
the rotation axis of the Kerr source. 

The specific angular momentum of the source J/M = a appears in equa- 
tions Q and as a^, so that these equations are invariant under the 



^ = -^GM^^^—^. (10) 
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transformation a —a. This implies that the sense of rotation of the black 
hole is of no consequence for the phenomena described in this section. This 
is consistent with the double-jet structure observed in astrophysics. It is 
therefore sufficient to treat only the motion of particles moving outward 
above the north pole of the black hole. 

It is important to mention here that the main results of this work are 
independent of the precise nature of the massive source. For the sake of 
concreteness, we have assumed that the source is a Kerr black hole with 
a < GM. Moreover, we will concentrate on a jet clump launched from 
ro ^ GM and moving along the positive z-axis. However, equations (jSJ- 
(|10j) depend upon a in terms of c? jr^ ^ which is very small for r > tq ^ a. 
In fact, our main results are insensitive to the ratio a/{GM). That is, our 
main conclusions remain the same even if the source is spherically symmetric 
(a = 0) and the jet clump is launched along a radial direction. 

The tidal equation in principle contains an infinite series in powers 
of the quantity Z/TZ, where TZ is the effective radius of curvature such that 
IZ^^ = \k\. For regions of interest in this work k ~ —2GM/r^ and hence 

Starting from r = ro corresponding to Z = at T = 0, we integrate equation 
© with a given initial Vq = ^(0) > 0. For a highly relativistic jet clump, the 
initial Lorentz factor Fq = (1 — V^^)~2 is such that Fq ^ 1. The integration 
time T must be limited such that Z[T) ~ T < IZq, where IZq is given by 
equation at r = rg. It is important to note that the endpoint of the 
integration is only known implicitly in terms of the Boyer-Lindquist radial 
coordinate r. There is, of course, a definite relationship between the Boyer- 
Lindquist coordinates and the corresponding Fermi coordinates. However, 
this relationship remains implicit here; making it explicit would immediately 
imply that we know the tidal equation Q to all orders 

The first-order approximation of the tidal equation corresponds to the 
generalized Jacobi equation, 

g + A;(l-2y2)z = 0. (12) 

This equation has solutions corresponding to V = ibl/-v/2, or F := (1 — 
y^)^2 = -^2, where the jet moves with uniform speed relative to the refer- 
ence particle. For our reference particle with 7 = 1, these special solutions 
are attractors |2] • Starting near the Kerr black hole and moving along 
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Figure 1: Plot of the Lorentz factor T versus T/{GM) based on the inte- 
gration of equation (|12j) with initial data ro/(GM) = 10, Z{Q) = and 
F(0) = y/rf^/ro for a/{GM) = 1. The graph illustrates the decelera- 
tion of clumps with initial Lorentz factors Tq = 1000, 100 and 10. 
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Figure 2: Plot of the Lorentz factor F versus T/{GM) based on the inte- 
gration of equation H12|) with initial data ro/{GM) = 100, Z{0) = and 
V{0) = v^Fg - 1/Fo for a/{GM) = 0. The graph illustrates the decelera- 
tion of clumps with initial Lorentz factors Fq = 1000, 100 and 10. 
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the positive z-axis, an initiahy ultrarelativistic jet clump decelerates rapidly 
relative to the reference particle (7 = 1) according to equation (|12() . This 
initial deceleration and approach to l/\/2 is depicted within the domain 
of validity of equation (|12() in figures ^ and HI for three clumps with initial 
Lorentz factors of Tq = 1000, 100 and 10. 

The critical speeds ±l/-v/2 that appear in equation ((T^ correspond to 
the critical Lorentz factor Tc = \/2. A particle of mass m with a Lorentz 
factor of Tc has kinetic momentum m, hence its de Broglie wavelength equals 
its Compton wavelength. For our present purposes, the motion of a particle 
is ultrarelativistic if its corresponding Lorentz factor F exceeds Tc- For 
F <C Tc, the tidal equation has already been extensively studied in the field 
of a Kerr black hole [31 [H] • 

The novel deceleration effect is thus a consequence of equation ((2)) re- 
stricted to ultrarelativistic motion along a certain straight line. It is there- 
fore interesting to consider the application of the generalized Jacobi equation 
to the motion of astrophysical jets ^I21- Thus imagine a rotating gravita- 
tional source with a double jet system moving outward along its axis of 
rotation. We have assumed that the Fermi coordinates established along an 
ambient trajectory that lies on the axis of rotational symmetry of the source 
are such that the axis of rotation is the Z-axis. Concentrating on an ultra- 
relativistic jet moving along the positive Z-axis, one finds that equation 
in this case involves a remarkably strong initial deceleration followed by a 
gradual decline toward a terminal speed given by l/\/2. Starting with an 
ultrarelativistic jet, this terminal speed is approached asymptotically after 
an infinite period of time. On the other hand, the generalized Jacobi equa- 
tion is only valid for a rather limited timescale that nevertheless covers the 
initial rapid drop in velocity. 

Figures ^ and |21 illustrate the main result of this paper regarding tidal 
deceleration; therefore, it is worthwhile to develop an approximate analytic 
understanding of how this tidal deceleration comes about. To fix our ideas, 
we focus attention in what follows on a jet clump with an initial Lorentz 
factor Fq ^ 1 moving along the positive z-direction relative to a free test 
particle of the ambient medium. At T = 0, Z = 0, i.e. the jet starts from 
the position of the reference particle on the z-axis, where 

(^r + ^^2GMr)'^' (13) 
according to equation Q with 7 = 1 and r(T = 0) = tq 3> 2GM. It follows 
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that in equation @, 

GM , , 
■ 2- (14) 

Starting from equation l|12|) and writing (fiZ/dT"^ = VdV/dZ, it is sim- 
ple to show that equation ((T^ can be completely integrated for a constant 
curvature k; in fact, the phase diagram for this case has been given schemat- 
ically in figure 2 of yLj. For k = k{T), however, let us note that 

d , /1-2V^\ 



In 



Since /c < 0, it follows form equation (|15j) that if Vq < 1/2, V'^ increases 
toward 1/2, while if Vq > 1/2, V'^ decreases toward 1/2. Integrating equa- 
tion (|15|) . we find 

In (j^^) = 4 £ k{T')V{T')Z{T')dT'. (16) 

To find an approximate solution of this equation we write equation IjlGj) in 
the form 

ln(^^-^J=4y^ k{T')V{T')Z{T')dT'. (17) 

Assuming that F » 1 and Fq ^ 1, the left side of equation ((T7|) can be 
approximately written as 2(Fq ^ — F~^). Next, in the integrand of equation 
(|T7|) we let V{T') k, 1 and Z{T') w T', and note that 



T 



k{T')T'dT' 



ln(l + A) + 



1 



1 + A 

using equation dJ. Here A is defined by 



(18) 



so that 2A/3 = T/TZq, using equation (|11|) . For A ^ 1, one can show that 

Hi + >^) + Y^^i + l>^' + o{x'). (20) 
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Equations (HTJ-dini) thus imply that 



r-2«ro2+(^)- (21) 



or V'^ Vq — {T/TZq)^. This result provides an approximate solution of 
equation (|12() for Fq ^ 1, F » 1 and T <^ 2TZo/3. Indeed, a graph of 
equation (|2H) is almost indistinguishable from the corresponding numerical 
solution of equation ((T^ presented in figures ^ and [2J 

To get an estimate of the time T that it would take for a highly relativistic 
jet clump with Fq 3> 1 to decelerate to F, where Fq ^ F^ S> 1, we note that 
A « 3/(2F) once we neglect Fq^ compared with F~^. Thus equation pTj) 
can be written as 

^ ~ i /I ( ^0 Y (22) 
GM ^ fV 2 \Gm) ' ^ ' 

It follows that the deceleration time is essentially independent of Fq. This 
remarkable result is already evident from the form of equation H17|) : for 
Fq ^ F 1, the left side of equation (fT7|) is nearly independent of Fq. In 
figure n for instance, tq = 10 GM and the time it takes for F to decrease to 
F 5 from any initial Fq 5 should be k. 2\/5 GM based on equation \12\ . 
in good agreement with the numerical results of figure ^ Thus no matter 
how relativistic the clump may be initially at tq = 10 GM, it decelerates to 
F « 5 in a short time of order GM. Moreover, this time is within the limit 
of validity of the generalized Jacobi equation, since equation may be 
written as 

1^ - i « 1. (23) 

The difference between figureEand figure|2lis a reflection of the fact that the 
curvature k decreases as the inverse cube of the radial distance (r » GM) 
from the black hole and, using equation (|22() for a fixed F, oc Tq, which 
is reminiscent of Kepler's third law. 

Thus far we have considered only one reference particle of the ambient 
medium with 7 = 1. Let us now imagine such particles all with 7 = 1 
distributed along the rotation axis of the source and constituting "fixed" fea- 
tures of the ambient medium that can be used for reference purposes. More 
specifically, let the Boyer-Lindquist radial coordinates r^, z = 1, 2, . . . , A^, of 
the particles be such that ri < r2 < . . . < rj < rj+i < . . . < r^v; moreover, 
each interval is small enough that the generalized Jacobi equation is valid 
within this interval to a reasonably good approximation. That is, N differ- 
ent Fermi frames can be established that are centered on the N particles 
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and the generalized Jacobi equation can be applied separately within each 
interval [rj, rj+i], so that if the initial speed of the jet clump relative to rj is 
more than the terminal speed, the clump continues to decelerate within this 
interval, and so on. This is a direct consequence of the remarkable feature 
of equations (j5|)- (|lUj) noted above: they do not explicitly depend on any 
specific reference particle with 7 = 1. In this sense, the clump decelerates 
with respect to the ambient medium toward the terminal speed. As the jet 
clump approaches the terminal speed, the first-order tidal term proportional 
to 1 — 2V^ in equation Q approaches zero and hence the higher-order tidal 
terms, neglected in equation ©, become dominant. A complete knowledge 
of these terms is then required to predict the subsequent evolution of the 
flow. It may be that at such a time the clump is sufficiently far from the 
black hole that tidal effects become unimportant and plasma forces take 
over the dynamics of the clump [2]. 

Consider now the tidal equation Q in the nonrelativistic limit \V\ <C 
1. The first-order approximation leads in this case to the standard Jacobi 
equation, i.e. equation ((T^ but with 1 — 2V^ replaced by unity. As explained 
before, for a clump along the rotation axis at r with ^ 3a^, k ~ —2GM/r^ 
just as in Newtonian mechanics. The tidal force of the central source in this 
case simply tends to pull the clump and the reference particle apart, leading 
to a relative tidal acceleration. The result of the integration of the standard 
Jacobi equation is in this case essentially the same as in the Newtonian 
theory 0101) which is studied in the next section. 



3 Comparison with Newtonian gravitation 

It is useful to compare the main results of the relativistic approach with the 
Newtonian theory of tides, since our main conclusions regarding tidal decel- 
eration and terminal speed have no counterparts in the Newtonian approach 
and hence violate our "nonrelativistic" intuition. 

Let us consider in the background inertial coordinate system {t,x,y,z), 
the Newtonian equations of motion of the jet clump at r and the reference 
particle at rp along the rotation axis of the source 

(fr _ _GM_ cPrp _ _GM_ 



We are interested in the relative motion; therefore, with = r 

(fQ GM , GM 

p 



dt^ {rp + Cr^ rl ■ ^^^^ 
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For C < '^pj tliis equation can be written as 



^ = ^C + 0(C), (26) 

which should be compared and contrasted with equation @. The higher- 
order Newtonian tidal terms are considered in [Appendix B 



Let us assume as before that the reference particle has zero energy, so 
that its kinetic and potential energies add up to zero; then, {drp/dt)"^ = 
2GM/rp can be simply integrated and the result is 

rl'^+^-V2GMt)"\ (27) 

where = ro at t = 0. Substituting this result in equation (|25|) . we find to 
linear order in C, 

(fC 2GM , , , 

^ rC = 0, (28) 



dt2 [3/2 



which can be solved with boundary conditions that C = dC,/dt = vq 
at t = 0. Here > is the initial speed of the jet clump relative to 
the reference particle. The radial Newtonian tidal force tends to pull the 
particles apart; therefore, ( = dC,/dt monotonically increases for the duration 
of the validity of the approximation used here. This is the complete opposite 
of what happens in the relativistic theory for an initially ultrarelativistic 
relative speed. 

An interesting feature of equation ()12() is therefore that the nature of the 
tidal force changes depending on whether is below or above the critical 
value = 1/2. In fact, for < 1/2 the tides behave essentially as in 
Newtonian mechanics, i.e. the particle accelerates — its speed approaching 
l/\/2 very slowly — relative to the reference particle. The character of the 
tides in a model star falling into a Kerr black hole has been investigated 
in detail in [31 |S] using the Jacobi equation; the results are similar as in 
Newtonian gravitation except possibly when the model star is very close to 
the horizon of the black hole. On the other hand, for V'^ > 1/2 the tides 
behave in the opposite way (see figures ^ and [JJ . 

Comparing the Newtonian behavior with the consequences of equations 
H12|) ~ ()14|) . we recognize the significance of the relativistic tidal accelera- 
tion term 2kV'^Z for the difference between the results of the two theo- 
ries. Thus tidal deceleration and terminal speed exhibited by ((T^ are due 
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to this gravitomagnetic (i.e. motional) tidal effect that is purely relativis- 
tic and thus goes beyond Newtonian gravitation. The higher-order tidal 
terms in equation © are treated using a certain post-Newtonian 
approach in Appendix B 



^-^0101 


= '^-^0202 = 


~2^' ^'^0303 


^-^2323 


= '^-Rsisi = 


-k, '^i?l212 = 


^-^0123 


= '^^0231 = 


~2^' ^-^0312 



4 Tidal dynamics 

It is interesting to explore the dynamics of the three-dimensional general- 
ized Jacobi equation ((2)) using the same reference trajectory © as before. 
The components of the curvature tensor in the Fermi coordinate system 
established along this worldline are given by |Hj 

k, (29) 
-k, (30) 

q; (31) 

these are the only nonzero components except, of course, for the symmetries 
of the Riemann tensor. Here (0, 1, 2, 3) = (T, X, Y,Z), k is given by equation 
(|l()j) and q is directly proportional to the angular momentum of the source 
and is given by 

37.2 _ „2 

q = 2GMa-^ ^. (32) 

(r2 + a2)3 ^ ^ 

One can interpret q as the "magnetic" curvature, while k has the interpreta- 
tion of "electric" curvature. It is a significant feature of equations H29() - (|5T|) 
that they are independent of the Lorentz factor of the reference trajectory 
7. This is a result of the fact that the axis of rotational symmetry in the 
Kerr metric provides two special tidal directions — corresponding to ingoing 
and outgoing trajectories — such that k and q are independent of any Lorentz 
boosts along these directions jHHSl^J. These special tidal directions are 
related to the repeated principal null directions of the curvature tensor and 
therefore the degenerate nature of Kerr spacetime. That is, the special tidal 
directions exist at each event in Kerr spacetime, since it is a vacuum space- 
time of type D in the Petrov classification jHl ll()| lllj . The significance of 
such directions in connection with the direction of an astrophysical jet has 
been pointed out in to preserve the collimation of a highly relativistic 
jet against tidal disruptions caused by the central source, the most natural 
jet direction would be a special tidal direction. Moreover, of all the special 
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tidal directions the Kerr rotation axis has the additional advantage of being 
in conformity with the symmetry of the configuration under consideration. 
Using equations (pn |) - (|5T|) . it is possible to express equation ^ as 

X - ^kX[l - 2X^ + ^(2y2 _ ^2^] + ^kX{5YY - 7ZZ) 

+q[XYZX - ZY{1 + X^) - 2YZ] = 0, (33) 

Y - ^kY[l - 2Y^ + ^(2^^ - Z^)] + ^kY{5XX - 7ZZ) 

-q[XYZY - ZX{1 + _ 2XZ] = 0, (34) 

Z + kZ[l- 2Z^ + -(X^ + y2)j + lkZ{XX + YY) 
3 3 

-q{XY - XY){1- Z"^) (35) 

where X = dX/dT, etc. These equations reduce to the one-dimensional 
generalized Jacobi equation (|12|1 for X = Y = 0. It is simple to check that 
equations H33|) - H35() are invariant under the transformations {X, Y, Z; k, q) i— > 
(y, X, Z; k, —q) and {X, Y, Z; k, q) i— > {X, Y, —Z; k, —q), since q is directly 
related to the rotation of the source. Furthermore, equations (|33|) - H35() are 
also invariant under parity or time reversal if in each case we also change q 
to —q. 

It is interesting to note an important symmetry of the system (|33)1- 
(|35j) under rotations about the Z-axis. The rotational symmetry about 
the direction of the motion of the reference particle implies that there is a 
degeneracy in the choice of the X and Y axes of the tetrad frame of the 
reference observer. Once an orthonormal tetrad is chosen, however, it is 
then parallel transported along the reference worldline. Nevertheless, the 
axial symmetry of the configuration under consideration implies that under 
a rotation about the Z-axis with a constant azimuthal angle (j)o, 

X' = X cos (po + Y sin cpQ, (36) 
y = -Xsin(/.o + ycos0o (37) 

and Z' = Z, the system (|33 |) ~ (|35j) remains invariant. This symmetry can, of 
course, be checked directly; for instance, for (f)Q = tt/2 we have X' = Y and 
Y' = —X and it is simple to see that the transformation (X, Y) (Y, —X) 
leaves the system invariant. 

To express k and q as functions of the time T, we need an explicit 
solution for the trajectory (jSJ. To this end, we set 7 = 1 as before and 
integrate dr/dT = V 2GMr / \/r^ + a? with the initial condition that r = ro 
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at T = 0, where Tq > 3a^. We have thus far focused our attention on an 
ultrarelativistic jet clump that is launched from (X, Y, Z) = at T = 
exactly along the rotation axis of the source such that X and Y remain zero 
throughout the flow. What happens if the jet is not launched exactly along 
the Z-axis? To answer this question, we must study the linearization of 
equations (|5n |) - H35|) about a solution Z = Z{T) of equation l(T^ . 

To obtain the linearization of system (|33|) - (|35jl around the motion of an 
ultrarelativistic clump along the axis of symmetry, we keep only terms linear 
in X, X,Y and Y and find that 

X - ^kiZZ)X - h{l - \z'^)X = q{ZY + 2ZY), (38) 
Y - h{ZZ)Y - h{l - '^Z^)Y = q{ZX + 2ZX), (39) 

while the third equation simply reduces to equation 1)121) . 

At any given time T, the new system (jnH|) - (|39|) represents two damped 
harmonic oscillators that are coupled through the magnetic curvature g, 
i.e. the coupling is due to the rotation of the source; moreover, the system 
decouples if it is transformed to the new variables X ^Y . It is possible to 
obtain a simple analytic solution of this system that is linearized about the 
solution of equation dJ, which amounts to Z{T) k, TZow{1 — uP' /&) 
with w = T/1Zq. Writing 



1 



-:^[l-3w + 0iw% q ^ 3—-^[l - 4w + Oiw"^)] (40) 



for < U) < 1, we find from equations (|HH|) ~ (|39j) that 

X{T) = X(0)TF{^) + ^^Y{0)TH{^), (41) 
TZo 2 ro 7^o 

Y{T) = Y{Q)TF{^) + ^^X{0)TH{^), (42) 

where only terms linear in a/ro <^ 1 are taken into account and F and H 
are given by 

F{w) = 1- ^w'^ + lw'^ + 0{w^), H{w) = w^ -2w'^ + 0{w'^). (43) 
12 8 

These results, which are approximately valid for To 1 and T <^ TZq, 
provide insight into the general behavior of the flow near the jet that moves 
along the Z-axis. 
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Figure 3: Plot of the Lorentz factor F versus T / (GM) based on integration 
of equation ^ with initial data X = at T = with X{0) = 2^/6/5, 
3\/TT/10, -v/399/20 and 7\/5T/50 corresponding respectively to Fq = 5, 10, 
20 and 50. In this plot a/{GM) = 1 and ro/(GM) = 10. The graph 
illustrates acceleration of the particle such that F essentially approaches 
infinity at T/{GM) ^ 3.9, 1.8 and 0.7 for Fq = 10, 20 and 50, respectively. 
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Figure 4: Plot of the Lorentz factor T versus T/{GM) based on integration 
of equation (EH) with initial data X = at T = with X{0) = 2\/6/5, 
3\/TT/10, \/399/20 and 7-v/5T/50 corresponding respectively to Fq = 5, 10, 
20 and 50. In this plot a/{GM) = and ro/{GM) = 100. The graph 
illustrates acceleration of the particle such that F essentially approaches 
infinity at T/ (GM) w 55 and 21 for Fq = 20 and 50, respectively. 
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The axial symmetry of the system p3() - ()35() under rotation about the 
Z-axis impUes that it would be advantageous to express this system in cylin- 
drical coordinates. This is done in Appendix A} where an approximate con- 



stant of the motion for the system H33|) - ()35() is also derived based on axial 
symmetry. 

An important feature of equations (|33|) - H35() is that for Y = Z = 0, this 
system reduces to 

X - ^kX{l - 2X^) = 0. (44) 

This equation has the same form as equation (|12() except that the nega- 
tive curvature k in (|12|) is replaced by —k/2. It follows that a particle 
launched from {X, Y,Z) = along the AT-direction with a speed above 
the critical speed l/\/2 is accelerated as in figures 01 and [l] It is a con- 
sequence of the axial symmetry of the system (|33j) ~ (|35|) that the form of 
equation (|44() is generic for motion along any direction perpendicular to the 
Z-axis (see Appendix A ). 



Equation (|44() again contains the critical speed l/\/2. In this approx- 
imation, if the initial speed is below this critical value, then the particle 
starting with a positive initial velocity decelerates over time intervals that 
are much longer than those corresponding to the validity of the approxima- 
tion. While the asymptotic behavior of the system has no physical meaning, 
one can show that the trajectories in the phase plane spiral around the 
origin. 

If the initial speed is above the critical value, then the particle acceler- 
ates and its Lorentz factor approaches infinity in finite time. In the latter 
case, it is interesting to extend our analytic approximation scheme involving 
equations in section [3 for Tq > 1, T > 1 and T <C 27^o/3 to 

equation H44|) by the simple substitution k ^ —k/2. Instead of H21() . we get 



which implies that T diverges at 

Td 1 / ro \3/2 



GM To\GM 



(45) 



(46) 



As expected, is inversely proportional to Tq for a fixed tq. For Tq = 
10, 20 and 50 in figure equation ^ predicts that Td/{GM) k, 3.16, 
1.58 and 0.63, which compare favorably with 3.9, 1.8 and 0.7 obtained by 
the numerical integration used to produce the figure. Similar results are 
obtained for figure for instance, Td/{GM) ~ 50 and 20 for Fq = 20 
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and 50, respectively, which should be compared with 55 and 21 based on 
numerical results presented in figure [l] The approximation ()46() improves 
as Tq increases. 

The phase diagram for equation 1)441) in the case of constant curvature 
is given in figure 3 of PP. It follows from this equation that 

which can be obtained from (|16|) by i— > —k/2. The left side of equation (|47|) 
is positive; therefore, if Vq < 1/2, then V'^ decreases toward zero. Moreover, 
if V^^ > 1/2, then V'^ monotonically increases with time. 

It is important to compare and contrast figures |21 and |11 the curvature is 
weaker for larger tq; therefore, for a given initial Fq, it takes longer for the 
particle's speed to reach unity (F oo). In fact, Tj oc rg from equation 1)46(1 . 
as can be verified based on the data given for figures 01 and EI 



5 Tidal acceleration 

We have demonstrated in the previous section that tidal acceleration can 
occur in directions perpendicular to the Z-axis if the speed of the clump 
exceeds the critical speed of l/\/2. In the three-dimensional case, motion 
relative to the standard reference trajectory can occur in X, Y and Z direc- 
tions. The main factors responsible for reversing the Newtonian behavior of 
the tides along these directions are 

ux=[l-2X^ + ^{2Y^-Z^)]^^^, (48) 
UY=[l-2Y^ + ^{2X' -Z^)]^^^, (49) 
uz=[l-2Z' + ^iX' + Y')]^^^, (50) 

based on equations (|3ci|l - (|35|) . For instance, ux < would imply that the 
character of the motion along the X-direction could be opposite of that 
expected from the Newtonian gravitation theory. 

Observe that I'x and i^y cannot both be negative. This follows from the 
inequality ux + > 0, since 

^{ux + uy) = -X^- Y^)t=o + - Z^)t=o, (51) 



20 



z 




Figure 5: Schematic representation of the boundary cone separating the re- 
gions of relative deceleration and acceleration. Particles whose initial veloci- 
ties are inside the critical velocity cone decelerate; particles whose velocities 
are outside the cone accelerate. The critical angle Q is such that cos^ ^ — \ 
and hence P2(cos0) = 0, where P2{x) = (3x^ — l)/2 is the Legendre poly- 
nomial of the second order. This angle is known in physics as the "magic 
angle" and occurs in several different contexts such as, for instance, NMR 
spectroscopy. 
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where the right side is a sum of two manifestly positive terms. Moreover, 
inspection of equations (|1S |) " H5U|) reveals that 

ux + ^>0, UY + ^>0, i/z + 1 > 0. (52) 

On the other hand, it is in general possible to have vz < and either I'x or 
i^Y also negative. An example of this situation is provided by the motion of 
a particle starting from the origin at T = with velocity components given 
by X = v/19/60, y = and Z = V5/3. While the dynamics of this case IS 
complicated, numerical experiments show that T slowly decreases with time. 

It follows from these results that highly relativistic particles emitted 
near a gravitationally collapsed configuration are accelerated along direc- 
tions normal to the rotation axis and hence immediately leave the black 
hole environment, while highly relativistic particles emitted along the rota- 
tion axis are decelerated and appear as jets. If the collapsed configuration 
is nonrotating, then deceleration occurs along the radial direction and ac- 
celeration occurs in directions normal to the radial direction. 

We have performed a series of numerical experiments to delineate the 
regions of acceleration and deceleration of particles relative to the refer- 
ence particle. The results are represented in figure El For ro/(GM) = 100, 
numerical experiments show that 9 is independent of Tq and is given approx- 
imately by tan^ = ^/2, so that 9 ~ 54.7°. In our numerical experiments, 9 
remained approximately constant for Fq > 5 and 5 < ro/{GM) < 100. 

To see how this comes about from an approximate analytical point of 
view, we consider equations (|.S.3j) - p5|) in cylindrical coordinates, namely, 
equations ()A2|) - (|A4|) . It follows from our numerical experiments that the 
influence of the rotation of the source on the relative motion is generally 
negligible. To simplify matters, the rotation of the source can therefore be 
neglected and we then confine our considerations to the {p, Z)-halfplane, 
since 99 = is a solution of equation (|A3|) with q = 0. Thus the equations 
of motion reduce to 

p-^kp = -kQp, (53) 
Z + kZ = -kQz, (54) 

where —k>0 and 

Qp = Pip"" + Iz^) - \pZZ, (55) 
Qz = Z(^-(? - 2Z2) + ^ppZ. (56) 
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If we set Qp = Qz = 0, equations and reduce to the standard 
Jacobi equation; therefore, the special deceleration/acceleration phenomena 
discussed in this work come about as a consequence of the existence of 
nonzero Qp and Qz- 

Let us now consider a particle starting out at time T = from the 
position of the reference particle at p = Z = 0. We will predict the motion 
of the test particle, determined by equations and (jSlJ, in case its initial 
velocity vector is {p,Z) = (Vq sin0, Vq cos 6*), where < Vq < 1 and the 
initial velocity vector makes an angle 9 with respect to the Z-axis. To do 
this, we let the function a be defined by 

piT) = a{T)Z{T), (57) 

substitute equation ([57j) in equation ([53|l and use equation (|5l|) to obtain 
the differential equation 

3 1 

aZ + 2dZ - -akZ + -akZ'^{2ap - 7Z) = 0. (58) 
2 3 

The Taylor expansions 

a{T) = {l + aiT + ^a2T'^ + ---)tan9, (59) 
Z{T) = {Vo cos e)T + ■■■ , (60) 

M^) = -7^o-2(l- At + ...), (6i) 

/-Co 

may be employed in equation H58|) to show that 

ai = 0, 02 = -^^0 ^- (62) 

Thus for T <C TZq, neglecting terms of second order and higher in T/TZq 
compared to unity, the particle trajectory is such that 

p=Ztan6, p=Ztan9 (63) 

are approximately valid. Plugging these relations in equations H55|) and 1)56(1 . 
we find that for T <^ TZq, 

Qp = pZ^itan'^e -2), (64) 
Qz = ZZ^{tan^e -2). (65) 

For tan^ 6 < 2, we have deceleration, while for tan'^ 6 > 2, we have accel- 
eration according to equations (|53() and Thus the boundary region is 
given by tan^ 9 = 2, as illustrated in figure El 
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6 Discussion 



In this paper, we are interested in some of the observable consequences 
of general relativity involving ultrarelativistic flows near a gravitationally 
collapsed configuration. The measurement problem in general relativity is 
rather subtle. Suppose, for instance, that in the background of Kerr space- 
time one succeeds in solving certain covariant equations of magnetohydrody- 
namics (MHD) in Boyer-Lindquist coordinates. Though possibly interesting, 
this result per se may have little to do with understanding the observable 
consequences of general relativistic MHD in this case. As a first step, we 
therefore concentrate on an invariant characterization of the dynamics of 
free particles in the exterior Kerr spacetime; in particular, we study relative 
motion in the simplest possible situation, namely, the motion of free parti- 
cles relative to a certain fiducial class of particles moving freely on escape 
trajectories along the rotation axis of the Kerr source. The relative motion 
of two nearby free particles is then determined by the generalized Jacobi 
equation. 

The generalized Jacobi equation has been known for three decades [3], 
but its physical consequences for relativistic relative motion have only re- 
cently received attention An important issue regarding the generalized 
Jacobi equation is that for ultrarelativistic relative motion, the duration of 
validity of this equation is rather short. However, we have shown that even 
within this limit, for ultrarelativistic outward motion along the rotation axis 
of a Kerr black hole, there is a remarkably strong initial tidal deceleration 
relative to the ambient medium regardless of the value of Tq ^ 1. To go 
beyond this initial deceleration, let us note that a significant feature of our 
basic equations (|8|)- (fT(H) is that they do not explicitly depend on the refer- 
ence trajectory corresponding to a "fixed" marker in the ambient medium. 
One can therefore imagine a number of such markers along the path of the 
jet such that from one marker to the next, the integration of equation (|12|) 
describes the force-free jet deceleration to a good approximation. In this 
way, one may approach the terminal speed of the jet by referring the jet 
motion to different markers. So long as the jet speed relative to any nearby 
ambient marker is significantly greater than l/\/2, it tends to decrease to- 
ward this terminal speed. 

It is possible to view the terminal speed in another way. Along a radial 
direction far from a massive source, its tidal influence on two neighboring 
test particles is a tidal acceleration if the relative speed of the particles is sig- 
nificantly below I/V2. However, this influence turns to a tidal deceleration 
if the relative speed is significantly above 1 / \/2- 
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We have neglected plasma effects in this paper; in fact, our treatment is 
valid for a force-free plasma. If the plasma is not force-free, the inclusion of 
its effects would enormously complicate matters jH^]. In any case, tidal 
deceleration is most effective very close to the central source, where the jet 
originates. Moreover, it is expected to be more significant for microquasars, 
since tidal effects near a black hole generally decrease as the inverse square 
of the mass of the source. Thus the results of this paper may well play a 
significant role in the dynamics of jets, especially in microquasars Hj. More 
generally, this work may be of interest in connection with the astrophysics 
of energetic particles such as the ultrahigh energy cosmic rays. 

Let us now turn to the motion of a particle normal to the rotation axis of 
the black hole. Relative to the reference particle, the particle accelerates if 
its initial speed exceeds the critical speed 1 / \/2- For ultrarelativistic relative 
motion, this tidal acceleration can lead to ultrahigh energy particles accord- 
ing to the generalized Jacobi equation; more generally, particle acceleration 
occurs outside a cone of angle 6, where the polar angle 9 is measured from 
the Z-axis and tan^ ~ \/2. As the particle accelerates beyond the Fermi 
system, one can in principle use other overlapping coordinate systems to 
describe the motion over an extended period of time. On the other hand, 
the tidal effects of the black hole diminish rapidly with increasing distance 
and hence the influence of the black hole can be neglected once the particle 
has moved sufficiently far away. 

The rotation of the source leads to a preferred radial direction, i.e. the 
axis of rotation, but is otherwise of little importance for the tidal dynam- 
ics presented in this paper. This remark applies equally well to the exact 
nature of the source, except that our results become physically significant 
near gravitationally collapsed objects, where tidal accelerations can be sub- 
stantial. 

The nature of the accretion mechanisms that create ultrarelativistic par- 
ticles near black holes is beyond the scope of our investigation. The analysis 
presented in this paper based on the generalized Jacobi equation indicates 
that once such particles are created near the poles of the black hole, those 
propagating mainly along the rotation axis of the black hole decelerate with 
respect to the ambient medium, while those propagating mainly normal to 
this axis can accelerate to almost the speed of light. 

The loss of kinetic energy for the motion of an initially ultrarelativistic 
particle along the rotation axis and the gain in kinetic energy for motion nor- 
mal to this axis may be attributed, in analogy with Newtonian gravitation, 
to the change in the gravitational potential energy, which is the tidal energy 
in the case under consideration in this paper. The underlying assumption 
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here is that the perturbation of the background geometry caused by the 
motion of the test particle may be neglected. This breaks down, however, 
when the particle is accelerated to almost the speed of light; therefore, the 
back-reaction on the motion is expected to moderate this singularity lead- 
ing to an ultrahigh energy particle. More generally, the gravitational and 
electromagnetic (in case of charged particles) radiations emitted by the de- 
celerating/accelerating particles should be taken into account as well. The 
tidal acceleration mechanism may provide the key to the explanation of 
recent Chandra X-ray observations of the Crab Nebula. 

It is well known that primary ultrahigh energy cosmic ray protons with 
energies above 10^*^ eV from extragalactic sources are expected to inter- 
act with the cosmic microwave background photons resulting in photopion 
production and pair creation j^^^. Therefore, it may be reasonable to as- 
sume that the most energetic particles that reach the solar system originate 
within our galaxy (see p3| for a recent review of cosmic ray astrophysics). 
In view of the foregoing results, it would be interesting to determine whether 
the directional distribution of the ultrahigh energy cosmic rays is correlated 
with the known microquasars, taking due account of the directionality of 
the acceleration zones normal to the jet directions. 



Appendix A Cylindrical coordinates 

In system (|5S|) - ()35|) . let us introduce the cylindrical coordinates p, ip and Z, 
where 

X = pcosip, Y = psinip. (Al) 
The resulting system is given by 

p-pif^- h[p{l - 2f? + ^p^^^ - + ^pZZ] 

+qp^{ppZ - 2Z) = 0, (A2) 
p(p + 2p^ + \kp^{bpp - 7ZZ) + q[2pZ + pZ{l + pV^)] = 0, (A3) 

Z + k{Z[l + i(p2 + /^V') - 2^2] + '^ppZ} 

-gpV(l - Z^) = 0. (A4) 

This system depends on ip only through 93 and (p and is therefore invariant 
under 1— > 99 + constant, which again reflects the axial symmetry of the 
system. 
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A remarkable feature of system HA2|) - (|A4|) is that for ip = and Z = 0, 
the system reduces to purely radial motion given by 

p - ^kp{l - 2p2) = 0, (A5) 

which is the same form as equation (|12)) except for k i-^ —k/2. This rever- 
sal of sign of the curvature leads to the acceleration of an ultrarelativistic 
particle that is launched from p = in a direction normal to the Z-axis as 
illustrated in figures 01 and ^ 

It has been shown in [jLj that the tidal equation is derived from a La- 
grangian given in Fermi coordinates by 



L = -^J-goo - 2goiV' - gijV'V^ , (A6) 

where the metric components in Fermi coordinates are given by equation 
Using equations ()29() - H31|) . it is possible to show that in our approximation 
the Lagrangian reduces to 

L = -^1 - X2 - y2 _ ^2 _ ^k{T)S - 2q{T){XY - XY)Z , (A7) 
where S{X, Y,Z,X,Y,Z) is given by 

5 = x^ + Y^-2Z^ + -{XY-XYf--{XZ-XZf--{YZ-YZf. (A8) 
3 3 3 

Expressed in cylindrical coordinates, L takes the form 



P^ = — (All) 



L = -y 1 - p2 _ p2^2 _ ^2 _ ifc(T)^ - 2q{T)p'^^Z , (A9) 

where 

^ = p2 _ 2^2 ^ 1^2(2^2 _ ^2)^2 _ _ pZ)\ (AlO) 

Here (/? is a cyclic coordinate; therefore, 

dL 

dip 

is an approximate constant of the motion for system (|A2|) - HA4|) given by 

p ^ [l + \k{T){2p^ - Z^)]p^^ + q{T)p^Z ^^^^^ 
^ ^1 - p2 _ p2^2 _ ^2 _ 1 A;(r)^ - 2q{T)p^pZ 

One can interpret P^p as a generalized orbital angular momentum about the 
Z-axis; its approximate constancy in time is related to the axial symmetry of 
the system (|A2|) - ()A4|) . We do not know if system ()A2|) - (|A4|) is Lagrangian. 



27 



25 



50 



75 



100 



125 



150 



175 



200 



Figure 6: Plot of the Lorentz factor T versus T/{GM) based on the inte- 
gration of equation HB1|) with initial data rQ/(GM) = 100, Z{0) = and 
V{0) = v^rg - l/To for a/{GM) = 0. The graph illustrates the decelera- 
tion of clumps with initial Lorentz factors Tq = 1000, 100 and 10. 

Appendix B Higher-order tidal accelerations 

It is in principle possible to obtain the infinite set of higher-order tidal 
accelerations that have been ignored in our analysis thus far. Moreover, we 
must impose the requirement that for a timelike geodesic (700 + '^dOiV^ + 
Qjky'^y^ < 0. The generalized Jacobi equation contains the first-order tidal 
terms [3]. Tidal terms of second and third order were first calculated in |S] 
and 0, respectively. Based on these previous detailed calculations, the 
first two such contributions to equation Q can be computed and hence the 
resulting augmented form of equation Q is 



+ \k'Z + -k"Z'^){l - 2V^)Z + -{K + -K' 
2 6 2 3 




(Bl) 
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Figure 7: Plot of the Lorentz factor T versus T/[GM) based on integration 
of equation (|BT3|) with initial data X = at T = with X(0) = 2\/6/5, 
3\/TT/10, \/399/20 and 7-v/5T/50 corresponding respectively to Fq = 5, 10, 
20 and 50. In this plot a/{GM) = and ro/{GM) = 100. The graph 
illustrates acceleration of the particle such that F essentially approaches 
infinity at T/ (GM) w 63 and 21 for Fq = 20 and 50, respectively. 
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Here V = dZ/dT, k is given by equation (jlUj) as before and 

k' (T) = ^-Ro303 ; 3 = RfMupa ; S A^gj AJ^g) A^^^ XJ^^^ A^g) , (B2) 

k" (T) = '^i?0303 ; 33 = Rpupa ; 5uj Aj^Q) A^^g) A^g^ A^g) A^g) A^'g) , (B3) 

K{T) = '^iio303 ; = Rpupa ; S X'^q^ A('g) A^q^ A^'g) A^g) , (B4) 

K'(r) = '^i?0303;30 = Rpvpa ;5ui A('o) AJ^s) A^q^ A[g) A(g) A^'g) , (B5) 

along the reference worldline where dx'^/dX^ = principle, the 

projections of the covariant derivatives of the Riemann curvature tensor of 
Kerr spacetime on the tetrad frame of the reference particle evaluated along 
its worldline are required for the determination of the higher-order tidal 
accelerations. However, for vq S> GM , it is simpler and more instructive 
to employ the linear approximation of general relativity as follows. The 
linearized Kerr metric in the post-Newtonian framework is given by 

{l-2^){dx' + dy^ + dz^), (B6) 



ds"^ = -{1 + 2^) dt^ + ^^{ydx - xdy)dt 



where <I> = —GM/r is the Newtonian potential. This metric can be written 
as Qpu = Vpu + h^y, where /iqo = -2^', hj = -2^6ij and 

h,, = -2G^^fi^ (B7) 

is proportional to the gravitomagnetic vector potential. The Riemann tensor 
in the linear approximation is given by 

Rpupa — 2 (^/'O" , >^P ^ hup , pa , pp ^PP , ua) i (-^^) 

which can be represented in the standard manner as a symmetric 6x6 matrix 
with indices that range over the set {01, 02, 03, 23, 31, 12}. For the linearized 
Kerr metric ()B6|) . we have such a representation in terms of symmetric and 
traceless 3x3 matrices E and B, 



E B 
B -E 



(B9) 



where E and B are the electric and magnetic parts of the curvature and 

GM x^x^ 

E,, = i?o.o, = = -^{S^, - 3^), (BIO) 
B^J = lejkeRoiM = -3^[xV^' + x^ j' + {6ij - 5^)r • J]. (Bll) 
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In the linear approximation under consideration here A'^j-^^ along the 
reference worldhne in effect reduces to 6^a in the computation of the Rie- 
mann tensor and its covariant derivatives (which reduce to partial deriva- 
tives) as measured by the reference particle. It is therefore interesting to 
note that along the reference trajectory {x = y = 0, z = r), the Riemann 
tensor measured by the reference particle turns out to be identical to equa- 
tions (|29|l - (|,'-{lj) without the terms, as can be verified by direct calculation 



using (p310|) and (|BTTJ). 

The rotation of the source plays a rather minor role in the numerical 
results of this paper; therefore, for the sake of simplicity we will set J = 
in the computation of k', k" , K and K' . We find that K = K' = {) and 
k = -2GM/r3, 

= ;c" = -24^. (B12) 

SO that k' = and k" = A;,33 along the reference path (z = r). It follows 
that the results in equation ()B12|) agree with the higher-order Newtonian 
tidal terms in equation H26() . Thus the first term on the right-hand side 
of I|B1|) contains the higher-order Newtonian tidal terms once 1 — 2V'^ i— > 1 in 
the nonrelativistic limit. Equation ()BH) can now be integrated with r = r(T) 
given by (jl^j) ; the results of this integration are presented in figure where 
the Lorentz factor versus T is given for ro/(GM) = 100 and initial Lorentz 
factors 1000, 100 and 10. This figure should be compared and contrasted 
with figure 121 

Let us now turn to the higher-order tidal corrections to equation 1)44(1 . 
Based on the detailed treatments in j^jQ, the result is 

= —X{1 — 2X^)(-^iioi01 + ^■^-^0101;!-'^ + ^■^-RoiOl;!!-'^^) 

+ 2"^^"'^(^-^oioi ;0 + 2 '^-^0101; 10-^) — 3"'^'^(-'^ + -^^)('^-Roioi)^ 
+ 0{X^). (B13) 

Using our approximation scheme with J = 0, we find that along the trajec- 
tory ^i?oioi ; = '^^oioi ; 10 = and 

Fr. _GM p -n - mi 4^ 

-KoiOl — — q-7 -Koioi;!— U, -Koioi;ll — — y — (,Di4j 

Equation HB13|) can be integrated with r = r(T) given by (|13|) : the Lorentz 
factor versus T is given in figure [71 for rQ/{GM) = 100. The results should 
be compared and contrasted with those given in figure 01 The terms in 
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equations (|B1|) and (|B13|1 proportional to k'^ make negligible contributions 
to the numerical work presented in figures Inland [3 respectively, in agreement 
with our approximation scheme. 

Within the limit of validity of the tidal equation, the higher-order tidal 
terms provide small corrections to the generalized Jacobi equation as in fig- 
uresiniandlTJ To demonstrate that the main results of this paper, namely, the 
strong initial tidal deceleration of ultrarelativistic flows along the rotation 
axis and the acceleration normal to this axis will remain unchanged by the 
addition of the higher-order tidal terms, we have integrated equations (jBll) 
and ()B13|) with the same initial data as in figures |21 and |1J respectively. The 
results are presented in figures El and [3 respectively, which appear to be 
almost identical to figures [21 and 0J respectively and demonstrate that the 
higher-order tidal terms slightly enhance the deceleration effect and tend to 
moderate the acceleration effect. For instance, in figure |21 the graph for ini- 
tial Tq = 100 reaches T ~ 4.63 at T = 100 GM, while in figure El it reaches 
r ~ 3.44. Nevertheless, figures (HI and [3 demonstrate that the inclusion of 
higher-order tidal terms will not change the main results of this work. 

Concerning the tails in figures (21 and we note that so long as the 
initial Tq is well above the critical value of \/2, a temporal interval of inte- 
gration can be found such that the fiow decelerates as described in detail in 
this paper. However, the size of this temporal interval of validity decreases 
as the fiow speed approaches the critical speed. Thus sufficiently close to 
the critical speed our approach loses its validity. Moreover, based on the 
available terms of the tidal equation ()BH) . it is clear that the critical speed 
l/\/2 is not a feature of the general tidal equation in this case. It follows 
that our treatment breaks down when the jet clump decelerates to a speed 
sufficiently close to this critical speed. On the astrophysical side, we note 
that by the time the clump decelerates close to the critical speed relative to 
the ambient medium, the black hole may be sufficiently far from the clump 
that tidal effects may no longer play a significant role in the subsequent 
motion of the clump. 
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